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' Abstract. We consider the Shaped Partition Problem of partitioning n given vec- 

tors in real /c-space into p parts so as to maximize an arbitrary objective function 
' which is convex on the sum of vectors in each part, subject to arbitrary constraints 

on the number of elements in each part. In addressing this problem, we study the 
Shaped Partition Polytope defined as the convex hull of solutions. The Shaped Par- 
\my , tition Problem captures A/'T^-hard problems such as the Max-Cut problem and the 

' Traveling Salesperson problem, and the Shaped Partition Polytope may have expo- 

(— I , nentially many vertices and facets, even when k or p are fixed. In contrast, we show 

that when both k and p are fixed, the number of vertices is polynomial in n, and all 
vertices can be enumerated and the optimization problem solved in strongly poly- 
nomial time. Explicitly, we show that any Shaped Partition Polytope has 0{n''^^^) 
vertices which can be enumerated in 0{n'^ ^ ) arithmetic operations, and that any 
^ . Shaped Partition Problem is solvable in 0{n''P ) arithmetic operations. 

(N 

■ 1- Introduction 
O ! 

. The Partition Problem concerns the partitioning of vectors ^4^, . . . , A" in /c-space 

^ I into p parts so as to maximize an objective function which is convex on the sum 

I of vectors in each part. Each vector represents k numerical attributes associ- 

' ated with the zth element of the set [n] = {1, . . . ,n} to be partitioned. Each or- 
dered partition tt = (tti , . . . , vTp) of [n] is then associated with the k x p matrix 



A" 



J2ieTTi ^\ ■ ■ ■ : SieTTp whose jth column represents the total attribute vector 



^ , of the jth part. The problem is to find an admissible partition vr which maximizes 

^ I an objective function / given by /(vr) = C{A^) where C is a real convex functional 

on IR^^^. Of particular interest is the Shaped Partition Problem, where the admis- 
sible partitions are those tt whose shape (|7ri|, . . . , |7rp|) lies in a prescribed set A of 
admissible shapes. In this article we concentrate on this later situation. 

The Shaped Partition Problem has applications in diverse fields that include cir- 
cuit layout, clustering, inventory, splitting, ranking, scheduhng and reliability, see 
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[3, 4, 6, 11] and references therein. As we demonstrate later on, the problem has ex- 
pressive power that captures ATP-hard problems such as the Max-Cut problem and 
the Traveling Salesperson problem, even when the number p of parts or attribute 
dimension k are fixed. 

Our first goal in this article is to show that a polynomial time algorithm for the 
Shaped Partition Problem does exist when both p and k are fixed. This result remains 
valid when the set A of admissible shapes and the function C are presented by oracles 
(see Section 4 for the precise statement) : 

• Theorem 4.6: Any Shaped Partition Problem is solvable in oracle time 0{n''P^) . 

Since C is convex, the Shaped Partition Problem can be embedded into the problem 

of maximizing C over the Shaped Partition Polytope defined to be the convex 
hull of all matrices A'^ corresponding to partitions of admissible shapes. The class 
of Shaped Partition Polytopes is very broad and generalizes and unifies classical 
permutation polytopes such as Birkhoff's polytope and the Permutohedron (see e.g. 
[15]). Its subclass of bounded shape partition polytopes with lower and upper bounds 
on the shapes was previously studied in [1], under the assumption that the vectors 
A^, . . . , A" are distinct. Therein a procedure for testing whether a given A" is a vertex 
of was obtained. This procedure is simplified and extended in [10] . A related but 
different important generalization of classical permutation polytopes, arizing when 
algebraic (representation-theoretical) constraints, rather than shape constraints, are 
enforced on the permutations involved, was studied in [15] and references therein. 

Since a Shaped Partition Polytope is defined as the convex hull of an implicitly 
presented set whose size is typically exponential in the input size even when both 
p and k are fixed, neither a presentation as the convex hull of vertices nor as the 

intersection of halfspaces is readily available for it. Our second objective is to bypass 
this difficulty by proving that, nevertheless, for fixed p and k, the number of vertices 
of Shaped Partition Polytopes is polynomially bounded in n, and it is possible to 
explicitly enumerate all vertices in polynomial time (see Section 4 for the precise 
statements) : 

• Theorem 4.3: Any Shaped Partition Polytope has 0{n''''~^^^^^^) vertices. 

While this bound suffices for our algorithmic considerations in this article, we point 

out that in [9] we establish the more accurate bound of 0{n''^^^) on the number of 
vertices. 

• Theorem 4.5: All vertices of can be enumerated in oracle time 0(n'^^^^). 

In particular, it is possible to compute the number of vertices efficiently. This might 
be extendible to the situation of variable k and p, while counting vertices is generally 
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a hard task (sec [12, 13]), as is counting partitions with various prescribed proper- 
ties (see [2]). The vertex counting problem for variable k and p will be addressed 
elsewhere. 

The article is organized as follows. In the next section we formally define the 
Shaped Partition Problem and Shaped Partition Polytope. We demonstrate that the 
Max-Cut problem and Traveling Salesperson problem can be modeled as Shaped Par- 
tition Problems with fixed p = 2 and k=l, respectively, and that the corresponding 
polytopes have exponentially many vertices. In Section 3 we establish a sequence of 
lemmas in preparation for Section 4. We extend earlier results of [1] and provide a 
necessary condition for to be a vertex of V^. This leads to the study of parti- 
tions of points in space whose convex hulls are pairwisc separable (in the special case 
k = p = 2 such partitions had been studied quite extensively, see e.g. [14]; the case 
k — S,p — 2 has also been considered quite recently in [5]). In the final Section 4 we 
reward our preparatory results of Section 3 and estabhsh Theorems 4.3, 4.5, and 4.6. 

2. Shaped Partition Problems and Shaped Partition Polytopes 

For a matrix A, we use A^ to denote its ith column. A p-partition of [n] := 
{1, . . . , n} is an ordered collection tt = (tti, . . . , tt^) of p (possibly empty) sets whose 
disjoint union is [n]. A p-shape of n is a tuple A = (Ai . . . , Ap) of nonnegative integers 
Ai . . . , Ap satisfying Aj = n. The shape of a partition n is |7r| := (|7ri|, . . . , Itt^I). 
If A is a set of |>-shapes of n then a A-partition is any partition tt whose shape |7r| is 
a member of A. 

Let ^4 be a real k x n matrix. For each |>-partition tt of [n] we define the A-matrix 
of TT to be the k x p matrix 



A"" ^ 

where X^j^Tr := whenever ttj = 0. We consider the following algorithmic problem. 



Shaped Pcirtition Problem. Given positive integers p,k,n, matrix A e IR*^^", 
nonempty set A of p-shapes of n, and objective function on A-partitions given by 
/(tt) = C{A^) with C convex on H*^^, find a A-partition tt* which maximizes /, 
namely satisfies 

f{TT*) = max{/(7r) : |7r| E A}. 



The complexity of this problem depends on the presentation of A and C and on 
whether real arithmetic computations or Turing machine computations on rational 
data are assumed. Pleasantly, our algorithms work in strongly polynomial time and 
can cope with minimal information on A and C. So both A and C may be presented 
by oracles. All that is needed of an oracle for A is that on query shape A it answers 
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whether or not A G A. All that is needed of an oracle for C is that on query it 
answers with C{A'^). 

Since C is convex, the Shaped Partition Problem can be embedded into the problem 
of maximizing C over the convex hull of 74-matrices of feasible partitions, defined as 
follows. 

Shaped Partition Polytope. For a matrix A G IR'^^" and nonempty set A of 
/^-shapes of n we define the Shaped Partition Polytope to be the convex hull of all 
A-matrices of A-partitions 

:= conv {A'' : |7r| G A} C IR'^^^. 

We point out that is the image of the Shaped Partition Polytope P^, with / 
the n X n identity, under the projection X A ■ X . In [8] this is exploited, for 
the situation where the function C is linear and A = {A : / < A < -u} is a set of 
bounded shapes, to solve the corresponding Shaped Partition Problem for all n, p, k 
in polynomial time by linear programming over P^, and to efficiently test for a given 
partition tt whether A^ is a vertex of V^- 

Next we demonstrate that, even when either k or p are fixed, the Shaped Partition 
Problem may be jVP-hard, and the number of vertices of the Shaped Partition Poly- 
tope may be exponential. Therefore, polynomial time algorithms for optimization 
and vertex enumeration may (and, as we show, do) exist only when both k and p are 
fixed. 

Example 2.1. McLx-Cut Problem and Unit Cube. Fix p = 2. Let A be the set 

of all 2-shapes of n, let k — n and let A be the nxn identity matrix. The A-matrices 
of A-partitions in this case are precisely all (0, l)-valued n x 2 matrices with each 
row sum equals 1. The Shaped Partition Polytope has 2" vertices which stand 
in bijection with A-partitions, and is affinely equivalent to the n-dimensional unit 
cube by projection of matrices onto their first column. Now let G — {[n],E) be an 
undirected graph, and define an objective function fa as follows: for a 2-partition 
TT = (tti, 712) of [n] let /^(Tr) = |{e G -B : |e fl 7ri| = 1}| be the number of edges in the 
cut of G defined by tt. Since each A" is a distinct vertex of V^, there is a convex C 
on IR'^^^ such that /gItt) = C{A^) for all tt. The Shaped Partition Problem in this 
case is precisely the Max-Cut problem of finding a largest cut in a graph hence is 

Example 2.2. Travelling Salesperson Problem and Permutohedron. Fix 

k = 1. Let p = n, let A = consist of the single shape 1" = (1, . . . , 1), and 

let A = [1,2, ... ,n]. The matrices A^ in this case are simply all permutations of 
A. The Shaped Partition Polytope has n\ vertices which stand in bijection with 
A-partitions, and is the so-called Permutohedron. Now let D G IR"^" be a symmetric 
nonnegative matrix, where Dij represents the distance between locations i and j. 



A POLYTIME ALGORITHM FOR SHAPED PARTITION PROBLEMS AND POLYTOPES 5 



Define an objective function f^, as follows. For partition vr = ({zi}, . . . , let 
foiT^) = — 12j=i Di^^i-^-^ be minus the length of the Traveling Salesperson Tour on 
[n] given by vr (with in+i := ii). Since each A" is a distinct vertex of V^, there is a 
convex C on IR" such that /z)(7r) = for all tt. The Shaped Partition Problem 

in this case is precisely the Traveling Salesperson Problem of finding a shortest tour 
hence is jVP-hard. 

We finish this section with an example demonstrating how our results enable a 
polynomial time solution procedure which bypasses the need for checking the expo- 
nentially many feasible solutions present. 

Example 2.3. Splitting. The n assets of a company are to be split among its p 
owners as follows. For i = 1, . . . ,p, the ith owner prescribes a nonnegative vector 
Ai = (Aj 1, . . . , Ai^n) with Yl]=i A-ij = 1, whose entries represent the relative values 
of the various assets to this owner. A partition tt = (tti, . . . , tt^) is sought which 
splits the assets among the owners as evenly as possible, and maximizes the i^-norm 

{J2^=i I SjGTTi Ajl'^)^ of the vector whose ith entry J^jeir^ ^i.j ^^e total relative value 
of the assets handed to the ith owner by tt. This can be formulated as a Shaped 
Partition Problem with k = p, the p x n matrix A whose rows are the Ai, the set A 
of all p-shapes A of n with [^J < Aj < \^] for all i, and objective function /q(7r) — 

Cg{A^) where is the convex function on IR^'^^ defined by Cq{M) = Ef=i 
For fixed p, by Theorem 4.6 we can find an optimal partition in polynomial time 
while the number of partitions is exponential f2(p"n^). We note further 
that for g = 1 the functional Ci is linear on JRF^^ hence our results of [8] apply and 
yield a polynomial time solution even when p is variable. 

3. PREPARATION 

In this section we establish a sequence of lemmas in preparation for the final section. 

Two finite sets U, V of points in IR^ will be called separable if there is a vector h E IR^ 
such that h^u < Ii^v for all u E U and v E V with u ^ v. It is not hard to see that 
U, V are separable if and only if their convex hulls conv{U) and conv(y) are either 
disjoint or intersect in a single point which is a common vertex of both. In particular, 
if U and V are separable then \U r\V\ < 1. 

Let ^ be a given k x n matrix. For a subset S C [n] let A^ = {A'' : i E S} he 
the set of columns of A indexed by S (with multiple copies of columns identified). 
A p-partition tt = (tti, . . . , TTp) will be called A-separable if the sets A^'' and A^'' 
are separable for each pair 1 < r < s < p. So tt is A-separable if for each pair 
1 < r < s < p there is a vector hr^s G IR*^ such that h^^^A^ < hJ^A^ for all i E Tir 
and j E tt^ with A^ ^ AK We have the following lemma which, for a matrix A with 
distinct nonzero columns, specializes to a result of [1]. 
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Lemma 3.1. Let A be a matrix in IR'^^", let A be a nonempty set of p- shapes of n, 
and let t: be a A-partition. If is a vertex ofV^ then tt is an A-separable partition. 

Proof. The claim being obvious for p = 1, suppose that p >2. Let be a vertex 
of V^. Then there is a matrix C G R,'^^^ such that the hnear functional on IR*^^^ 
given by the inner product (C, X) = ^^=i C^l ^1 is uniquely maximized over 
at A^ . Pick any pair 1 < r < s < p, and let /i^.s = — C^ . Suppose there are i & TTr 
and j e TTs with A^ ^ A^ (otherwise Ai^'^ and A"^" are trivially separable). Let tt' be 
the A-partition obtained from vr by switching i and j, i.e. taking tt^ := tt,. U {j} \ {i}, 
< := 7r,U{i}\{i}, and tx[ := nt for all t r, s. Then (A^')^ = {A'^Y+A^-A' 7^ (^'^j'' 
hence A'^' ^ A"" . By choice of C we have (C, A'^') < (C, A'^) and so 

- A') = (C^ - C-fiA^ - A') = f^(C*)^((A'^)* - (A'^'Y) = (C, - A""') > 0. 

This proves that A'^'' and A"" are separable for each pair 1 < r < s < p, hence tt is 
A-separable. □ 

Let A be a A; X n matrix. A p-partition tt = (tti, . . . , tt^) will be called A-disjoint if 
conv(A'^'') and conv(A'^=) are disjoint for each pair 1 < r < s < p. Also, for A e IR*^^'* 
let A be the matrix obtained by appending a new row [1, . . . , n] to A. 

Lemma 3.2. Let A G IR^^" and let A be as above. For each A-separable partition tt 
there exists an A-disjoint partition n of the same shape |7r| = |7r| satisfying A" = A-. 

Proof. Let tt = (tti , . . . , Tip) be any A-separable p-partition. As any 1-partition is 
both separable and disjoint, wc may assume p > 2. Let tt be the lexicographically 
first partition for which {A-»} = {A""^} for i = l,...,p, where {A^} stands for 
the multiset of columns of A indexed by S. More precisely, tt is the partition such 
that for all u G IR'^ the following hold: if we put Ui :— {j G TTj : A^ = u} and 
Hi := {j G TTj : A^ = u}, then |f/j| = \Ui\ for alH = 1, . . . ,p, and i < j for all i E U_r 
and j G U_s with 1 < r < s < p. By definition n has the same shape as tt and satisfies 

Now consider any pair 1 < r < s < p. Since tt is A-separable, there is a vector 
hr,s G IR^ such that hJ.,A^ < h'^,,A^ for all i G vr^ and j G tt^ with A* 7^ A^ . Let /i,. 
be the (A; -|- l)-vector obtained by appending a last coordinate e to hr^s, with e > 
sufficiently small. For i G tt^ and j G vr^ consider the quantity 

- A) = hl,{A - A) + eU - ^)■, 

if A* 7^ A-' then, since e is infinitesimal, this quantity is positive; if A* = A^ then, by 
choice of TT, i < j and this quantity is again positive. So conv(A-'-) and conv(A-s) 
are disjoint. It follows that tt is A-disjoint and the proof is complete. □ 
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For V e M*^ denote by v the (A; + l)-vector obtained by appending a first coordinate 
1 to V. An orientation of a liyperplane if is a presentation of H as the zero set H = 
{x e K'^ : h^x = 0} of a linear form (with h G IR^'*'^). With an orientation of H, the 
closed and open halfspaces H- , if < below it are distinguishable from the ones H- , 
above it. The hyperplane H spanned by an ordered list , . . . ^v^] of k affinely 
independent points in IR'^ will be oriented as H = {x & IR*^ : det[^;^, . . . , v'^, a;] = 0}. 

A matrix A will be called generic if its columns are in affine general position. Let 
A e m'^^" be generic with n > k. With each triple (J, J^) where J is a fc-subset 
of [n] and ( J^, J^) is a 2-partition of /, we associate two ^-disjoint 2-partitions tt^ 
and TT^ as follows. Let [A^] denote the list of columns of A indexed by / and ordered 
as in A. Let H := aff([yl.^]) be the hyperplane spanned by A^ oriented by the ordered 
list [A^]. Let /< := {i e [n] : A' e H<] and /> := {i e [n] : A' e H>}. Define 
7r< := (/<UJ<,/>UJ>) and7r> := (/>U J>, /<U J<). Both 7r< and 7r> are A-disjoint 
since H can be slightly perturbed so as to strictly separate A^^ and A^^. 

Lemma 3.3. Let A e IR'^''" be generic. The set of A-disjoint 2-partitions is the 
set of all 2-partitions if n < k and is the set of partitions associated with triples 

(/, J'^, J^) as above if n > k. 

Proof. Un < k then the columns of A are affinely independent hence conv{A'^^) and 
conv(yl.'^2) are disjoint for any 2-partition tt. Suppose now n> k and let tt = (7ri,7r2) 
be any A-disjoint 2-partition. Since conv(A'^^) and cony{A'^^) are disjoint, a standard 
separation result assures the existence of a hyperplane spanned by columns of A such 
that A'^^ is contained in one of its closed halfspaces and A'^^ in the other. Since A is 
generic, this hyperplane H contains exactly k columns of A, so that H = afT(A^) for 
some /c-subset / of [n]. If the orientation aff([yl.-'^]) of H is such that A"^^ C H- and 
A^^ C H- then tt is the partition tt^ associated with the triple (7, / n tti,/ fl 772), 
whereas if A^'^ C H- and A'^'^ C H- then tt is the partition tt^ associated with the 
triple (/, Jn7r2,/n7ri). □ 

Let now A e IR'^^" be an arbitrary matrix. Define for each e > a perturbation 
A{e) e IR''^" of A as follows: for i = 1, . . . , n let the zth column of A{e) be A(e)* := 
A^ -\- eMl where M| := [i, i^, . . . , z'^]^ is the image of i on the moment curve in IR'^. 

Definition 3.4. Let A G IR*'^" and p > 1. A p-partition tt of [n] will be called 
A-generic if it is A(e)-disjoint for all sufficiently small e > 0. The set of A-generic 
p-partitions will be denoted by 11^. 

Lemma 3.5. Let A e R,'^^" and p > 1. For all sufficiently small e > 0, the per- 
turbation A{e) is generic and the set of A{e)- disjoint p -partitions contains the set of 
A-disjoint p-partitions and is independent of e and equals H^. 

Proof Regard the entries of A(e) as hnear forms in e. Recall that v denotes 
the {k -\- l)-vector obtained by appending a first coordinate 1 to v E M!^. Suppose 
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first that n > k. Then each determinant det[A(e)*", . . . , ^(e)**] is a polynomial of 
degree k in e whose leading coefficient is the nonzero Van der Monde determinant 
det[M^", . . . , Ml''], hence has finitely many zeroes. So for all sufficiently small e > 
all of these determinants are nonzero hence A{e) is generic. For n < k the same 
argument applies by restricting attention to the first n rows of A. 

Next consider any A-disjoint p-partition tt. The claim being obvious for p = 1, 
suppose that p >2. For each pair 1 < r < s < p, there is a hyperplane which strictly 
separates A^"" and A^". Clearly, for all sufficiently small e this hyperplane also strictly 
separates A(e)'^'' and A(ey% so tt is also A(e)-disjoint. 

Finally, we show that for all sufficiently small e > the set of A(e)-disjoint p- 
partitions is the same. We start with the case p = 2. If n < k then by Lemma 3.3 
this set simply consists of all 2-partitions. So assume that n > k. Let I = {ii, . . . ,ik} 
be any A;-subset of [n] with ii < • ■ ■ < i^. For each i ^ I, consider the following 
determinant, which is again a nonzero polynomial of degree k in e: 

D\e) := det[A(e)'^ , . . . , AieY" , 1(e)'] = ^ . 

j=0 

For all sufficiently small e > the sign of D^{e) equals the sign of the first nonzero 
coefficient among Dq,...,DI and is either positive or negative and independent of 
e. Denote this generic sign by sign(D*). Then = {i G [n] : sign(D*) < 0} 
and = {i G [n] : sign(D') > 0} are independent of e and therefore for each 
2-partition (J^, J^) of /, the partitions tt^, tt^ associated with the triple (/, J^, J^) 
are independent of e. So the set of partitions associated with triples (/, J"^, J^) is 
the same for all small e > and by Lemma 3.3 this is the set of 74(e)-disjoint 2- 
partitions. Now consider any p-partition n with p > 2 arbitrary. Consider any pair 
1 < r < s < p, and restrict attention to the submatrix [^'^'•UTsj Qf ^ consisting 
of the columns of A indexed by tt,. U tt^. The same determinantal argument as in 
the paragraph above shows that for all sufficiently small e the 2-partition {7rr,7Ts) of 
TTr U 7is is either [A(e)'^'''-''^'']-disjoint or not independently of e. So conv(A(e)'^'') and 
conv(v4(e)'^'') are either disjoint or not independently of e. Since this hold for each 
pair 1 < r < s < p, we see that tt is either /l(e)-disjoint or not independently of e as 
claimed. □ 

Let p > 2. With each hst [tt^'* = {'H'i^.t^I'^) : 1 < r < s < p] of (^2) 2-partitions of 
[n] associate a p-tuple n = (tti, . . . , TTp) of subsets of [n] as follows: for i = 1, . . . ,p 
put 

Since vr^ C 7r['* and tt^ C 7:2^ for all 1 < r < s < p, the TTj are pairwise disjoint. 
If also Uf^j^TTj = [n] then tt is a partition, and will be called the partition associated 
with the given list. If moreover each n^'^ is A-disjoint, then so is tt. 
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Lemma 3.6. For A e JR^^"' and p > 2, the set of A- disjoint p-partitions equals the 
set of p-partitions associated with lists of A-disjoint 2-partitions. 

Proof Let tt — (tti, . . . , Wp) be an A-disjoint p-partition. Consider any pair 1 < r < 
s < p. Since conv(A'^'') and conv(74'^^) are disjoint, there exists an oriented hyperplane 
Hr,s which contains no column of A and satisfies A'^'~ C H^^ and A'^'' C H^g. Let 
^r,s ^^r,s^ ^r,s^ A-disjoint 2-partition defined by vrp' ':= {i G [n] : A^ e i/<J 

and 772^ :— {i e [n] : A^ e H^,s}- Let tt' be the p-tuple associated with the n'^'^. 
Then the tt' are pairwise disjoint and for i = 1, ... ,p we have 

Since [n] — Uf^iTTi C UfLiTr-, it follows that n — tt' is the ;)-tuple associated with the 
constructed list of A-disjoint 2-partitions. □ 

Lemmas 3.5 and 3.6 imply our final preparatory lemma. 

Lemma 3.7. For A e IR''''" and p>2, the set of A -generic p-partitions equals 
the set of p-partitions associated with lists of A-generic 2-partitions. 

Proof. By Lemma 3.5 there is an e such that Il\ equals the set of 74(e)-disjoint 
2-partitions and 11^ equals the set of A(e)-disjoint p-partitions. The lemma follows 
by applying Lemma 3.6 to A{e). □ 



4. Vertex Enumeration and Optimization 

We now use the facts established in the previous section to provide an upper bound 
on the number of vertices of any Shaped Partition Polytope. We then proceed to give 
efficient algorithms for the enumeration of these vertices and for the solution of the 
Shaped Partition Problem. 

Lemma 4.1. For any A e IR^^" and positive integer p, the number of A-generic 
p-partitions satisfies |n^| = 0{n^^). 

Proof. If p = 1 the bound is trivial so assume p>2. By Lemma 3.5 and Lemma 3.3 
applied to the perturbation A{e) with e sufficiently small, ii n < k then T{\ consists 
of all 2" = 0{2^) = 0{n^) 2-partitions, whereas if n > A; then it contains at most 
two associated partitions per triple (/, J<, J^) with I a A;-subset of [n] and (J<, J>) 
a 2-partition of /. So |n^| is at most twice the number 2*^^^^ of such triples. By 

Lemma 3.7, the set contains at most one partition per each list of Q partitions 
from n^, hence, as claimed, 

|n^l<|nll® <(2'+'f?l)®-0(^'®)- ° 
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For A e IR'^^"' and nonempty set A of |)-shapes let 11^ be the subset of of A- 
generic A-partitions. For A e IR*^^" let A be obtained as before by appending a new 
row [1, . . . , n] to A. 

Lemma 4.2. Let A e IR'^^" and let A be a nonempty set of p- shapes of n. Any 
vertex ofV^ equals the A-matrix A^ of some partition tt in Yl\. 

Proof. Consider any vertex A^ of V^- By Lemma 3.1, the A-partition tt is A- 
separable. By Lemma 3.2, there is another A-partition tt such that A^ — A- and tt 
is 74-disjoint. By Lemma 3.5, tt e 11^ hence is in 11^. □ 

We can now conclude an upper bound on the number of vertices of any Shaped 
Partition Polytope, which suffices for our purposes here. In [9] we use a more delicate 
argument which is based on treating the multisets of columns of A in the original 

space, without hfting to (A;-|-l)-space, and derive the more accurate bound of 0{n^^^) 
on the number of vertices. 

Theorem 4.3. For any A e IR'^^", positive integer p and nonempty set A ofp-shapes 
of n, the number of vertices of the Shaped Partition Polytope is 0{n^''~^^^^^^). 

Proof Let A be as above. By Lemmas 4.1 and 4.2 above, 

|vert(7'^)| < |n^| < |n^| = Oin^^+'^i")). □ 

Note that this theorem stands in contrast with the fact that typically the number of 
A-partitions is exponential in n even for fixed p. 

We next describe an algorithm for the enumeration of H^. 

Lemma 4.4. There is an algorithm that, given positive integers p, k, n and a matrix 
A G IR'^^", enumerates the set 11^ of A- generic p-partitions using 0{p^ ■ n'^^^^^^^) real 
arithmetic operations. For rational A and fixed p and k this is a strongly polynomial 
time algorithm. 

Proof. If p = 1 then let 11^ := {([n.])} consists of the single p-partition ([n]). 
Suppose now p > 2. First, we describe how to produce the set 11^ of A-generic 
2-partitions. For each e > let A(e) E IR'^^" be a perturbation of A obtained, as 
before, by setting A(e)* := A(e)* + eM^, and let A(e)* be obtained, as before, by 
appending a first coordinate 1 to ^(e)*. If n < A; then, by Lemma 3.3, 11^ consists 
of all 2" = (9(2'^) = 0(n'^) 2-partitions. Suppose now n > k. Let I — {ii, . . . ,ik} be 
any A;-subset of [n] with ii < ■ ■ ■ < i^. For each i ^ / let 

j=0 

As explained in the proof of Lemma 3.5, D\e) is a nonzero polynomial of degree k in 
e whose sign -D'(e) is the same generic sign sign(D*) for all sufficiently small e > and 
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equals the sign of the first nonzero coefficient among D^, . . . , D\. To compute these 
coefficients, substitute e = 0, . . . , A; and compute each of the numerical determinants 
D*(0), . . . , D^{k). Then solve the system of linear equations I]j=o ^] — ^ ~ 

0, . . . , in the indeterminates Dq, . . . , D^. Since the defining matrix of this system 
is a Van dcr Monde matrix hence nonsingular, D^, . . . , Dl are uniquely determined 
and so is sign Each of the numerical determinants -D*(e) can be computed, 
and the system of equations solved, by Gaussian Elimination which takes 0{k^) 
arithmetic operations. So each sign(D*) is computable in 0{k^) operations. Therefore 
— {i e [n] : sign(£)*) < 0} and — {i e [n] : sign(D*) > 0} are computable in 
0{nk^) operations, as arc the two 2-partitions and tt^ associated with each triple 
(J, J'^,J^). As there are 2^(^^^ such triples, the set 11^ of A-generic 2-partitions can 

be computed in 2*^^^^nfc^ = 0{n''~^^) operations. 

To construct 11^, produce, using Oip^n) operations, from each list of partitions 
from n^, the associated p-tuple tt. Then test if it is a partition (i.e. if Uf^j^TTj = [n]). 
As there are at most {2^{^)^^^ such lists, the total number of arithmetic operations is 

bounded hy p^ni^J^i^^)^^ = 0{p^ ■n'^^^^'^^) which subsumes the work for constructing 
Il\ and is the claimed bound. Lemma 3.7 guarantees that each A-generic p-partition 
of shape A is included. The claim about rational data follows from the fact that 
Gaussian Elimination admits a strongly polynomial time algorithm (see e.g. [16]). □ 

We can now provide the algorithm for the vertex enumeration of any Shaped Par- 
tition Polytope. It suffices for the algorithm that the set of shapes A (whose size may 
be exponential in the rest of the data) is presented by an oracle that on query shape 
A answers whether or not A e A. 

Theorem 4.5. There is an algorithm that, given positive integers p, k, n, matrix A e 
IR'^^", and an oracle presentation of a nonempty set A of p-shapes of n, enumerates 
the set of vertices of the Shaped Partition Polytope V\ using 0{n^ ^ ) real arithmetic 
operations and oracle queries. For rational A and fixed p and k this is a strongly 
polynomial time oracle algorithm. 

Proof. Let A G be as before. Use the algorithm of Lemma 4.4 applied 

to A to construct the set of A-generic p-partitions 11^. For each partition n in 
this set, test if it is a A-partition by querying the A-oracle about its shape |7r|. In 
this way, produce the set 11^ of A-generic A-partitions. Then construct the set of 
k X p matrices U := {A'^ : tt e 11^} with multiple copies identified. This set U is 
contained in and by Lemma 4.2 contains the set of vertices of V^. So u E U 
will be a vertex precisely when it is not a convex combination of other elements of 
U. Testing it is a Linear Programming problem. To solve it, we shall appeal to 
Garatheodory's theorem and exhaustive search rather than invoke an LP-algorithm. 
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This allows to deal with real data and, for rational data and fixed k,p, obtain a 
strongly polynomial time procedure. So to test if -u G f/ is a vertex of do the 
following. Let d := dim(aff(?7 \ {u})) < kp + 1. For each of the (j^^^^^ subsets 
{ui, . . . , Ud} oi U \ {u}, check if it is an affine basis of aS{U \ {u}) by checking if 
dim('Ui, . . . , Ud) — d; if it is, obtain the unique expression u — HiUi of u as affine 
combination of the Ui. If /xi, . . . , /x^ > then u is in fact a convex combination of the 
Ui hence not a vertex of So ?/, is a vertex of if and only if for each affine basis 
we get some jii < 0. To compute d, to compute dim('Ui, . . . , m^), and to solve for the 
/ij, use Gaussian Elimination again. 

By Lemma 4.4 applied to A, the set LI^ can be constructed in 0(71^*^+^^^ ) opera- 
tions. The number of oracle queries is as guaranteed by Lemma 

4.1, and the number of matrices in U satisfies \U\ < |n^| < |n^| = 0{n'^''^^^(^>). As 
the effort of constructing is absorbed in the vertex testing procedure, we have 
that the total number of arithmetic operations and oracle queries is bounded by 
\U\(^^^~l^{kp + 1)^ = 0{n^^^^) as claimed. The claim for rational data follows from 
the strongly polynomial time computability of Gaussian Elimination. □ 

Finally, we provide an algorithmic solution of the Shaped Partition Problem we 
started with. It suffices that the convex functional C is presented by an oracle that 
on query with vr a A-partition answers with C{A^). The oracle for C will be called 
M -guaranteed if C{A'^) is guaranteed to be a rational number whose binary encoding 
is no longer than M bits for any A-partition tt. 

Theorem 4.6. There is an algorithm that, given positive integers p, k, n, matrix 
A e IR'^^", an oracle presentation of a nonempty set A of p-shapes of n, and an 

oracle presentation of a convex functional C on IR*^^^, solves the Shaped Partition 
Problem of finding a A-partition tt* which maximizes /(vr) = C{A'^) using Oin^"^ ) 
real arithmetic operations and oracle queries. For rational A, fixed p and k, and 
M-guaranteed oracle for C , the running time of this oracle algorithm is strongly 
polynomial in n, M and the bit size of A. 

Proof. Let A be as before. Use the algorithm of Lemma 4.4 applied to A to con- 
struct the set of partitions 11^, using • n^'^'^^^i^)'^^) arithmetic operations. Then 
filter out its subset of A-generic A-partitions by querying the A-oracle on each 

of the |n^| = 0(n('=+^)(2)) partitions in H^. Since C is convex, there will be an 
optimal A-partition vr for the Shaped Partition Problem for which A^ is a vertex of 
V^, hence, by Lemma 4.2, an optimal n E TL\. For each tt G 11^ compute A^ (using 
0{kn) operations) and query the C-oracle for the value C{A^). Pick tt* to be the one 
attaining maximum value. The number of operations and queries to the A-oracle and 
C-oracle is bounded, as claimed, by 0{kn ■ n^'^+^Ka)) + Qip"^ ■ n^^^+^Kaj+i) ^ o{n^P^). 
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For rational data, since the C-oracle is M-guaranteed, each of the comparisons of 
its answers needed to determine the optimal tt* is doable in time linear in M. Since 
is strongly polynomial time computable for fixed p, k, the claim about strongly 
polynomial time solution of the Shaped Partition Problem follows. □ 
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